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Abstract of EP 0315545 (A1) 

Electroencephalographic signals (EEG) collected by 
a set of electrodes (10) arranged at predetermined 
locations of the scalp are converted into digital form 
and transmitted to a processing unit comprising: 
means for optimal estimation of a cross-spectral 
matrix for the EEG signals collected; storage means 
In which are recorded information relating to a 
propagation model and allowing characterisation of 
at least one source vector defining relationships, 
functions of position parameters of the or each 
source, in the frequency domain, between a signal 
emitted by a bipolar source or each source and the 
signals available on the scalp, the propagation 
model being a model with three discontinuities 
comprising the brain, the skull and the scalp; and 
means for processing the estimated and noise- 
cleansed cross-spectral matrix in order, on the one 
hand, to identify the cerebral activity sources 
contributing to the EEG signals collected and 
determine their number and, on the other hand, 
estimate the coordinates of each source identified. 
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@ Disposltif de st§r6oloca]Isation spatiaie de sources d'activitds cdr6braies. 

@ Des signaux electroenceplialograpiiiques (EEQ) recueiHis 
par un ensemble d^electrodes (10) disposees en des emplace- 
ments predetermines du scalp sont convertis sous forme 
numSrique, et transmis a une unite de traitement coniprenant : 
des moyens d'estimation optimale d'une matrlce Interspectrale 
des signaux EEG recueillis ; des moyens de memorisation dans 
lesquels sont enregistrees des informations relatives a un 
modele de propagation et permettant de caract6rlser au moins 
un vecteur source definissant des relations, fonctions de 
paramdtres de position de la ou ciiaque source, dans le 
domaine frequentiel, entre un signal emis par une source 
dipolaire ou chaque source et ies signaux disponibles sur le 
scalp, le module de propagation etant un modele a trois 
discontlnuites comprenant le cerveau, le crane et le scalp ; et 
des moyens de traitement de la matrice interspectrale estimee 
d^barrassee du bruit pour, d'une part, identifier ies sources 
d'activltes cerebrales contribuant aux signaux EEG recueillis et 
determiner leur nombre et, d'autre part, estimer les coordon- 
nees de chaque source identifiee. 
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Description 

Dispositif de stereolocalisation spatiale de sources d'activit6s c6rebrales 

La presents invention a pour objet un dispositif de stereolocalisation spatiate de sources d'activites 
c§rebrales, c'est-^-dire un dispositif permettant de iocaliser dans le volume cerebral, par rapport a un systeme 
5 de coordonnees spatiales, une ou piusieurs sources a I'origine d'activit§s 6iectdques c^r^brales spScrfiques. 

Un domaine d'appllcation particulier de f'invention est la localisation spatiale d*une activity c§r6brale 
pathologique, par exemple un foyer epileptique. 

II est connu de realiser une cartograpiiie cerebrale a partir de signaux electroenoephalographiiques (signaux 
EEG) recueillis au moyen d'eiectrodes implantees en differents endroits du scalp d'un patient. La cartographle 
10 cerebrale consiste a etabllr des cartes illustrant la repartition spatiale de la puissance d'emission a differentes 
frequences par des nuances de gris ou des couleurs differentes, 

Un dispositif d'analyse spectrale permettant d'obtenir de bons resultats en cartographle cerebrate a deja 
ete decrit dans la demande de brevet franpais No 87 07 600. Dans cette den tande, 11 est notamment decrit un 
dispositif d'estimation d'une matrlce Interspectrale de signaux Ir.EG recueilli^^ par des Electrodes, par mise en 
IS oeuvre, sur les signaux EEG, d'un a|gonthme recursif cte blanchiment et decorr6lation spatiale avec 
optimisation du compromis bials-variance pour chacun des 6!6ments de ia matrlce. 

La cartographle cerebrale permet de foumir des Informations precieuses quant a Texistence de sources 
d'activites cerebrales et a leur emplacement approximatif . Toutefois, en particulier lorsqu*une exerese doit 
etre pratiquee, 11 est necessaire de Iocaliser avec precision ia source d'activlte cerebrale pathologique. A 
20 f'heure actuelle, cette localisation est realisee dans la majorite des cas par implantation d*electrodes au sein 
de ia masse cerebrale, a travers ia borte cranlenne, ce qui est particullerement traumatisant. 

Aussi, la presente invention a-t-elle pour but de foumir un dispositif permettant, a partir d'une matrice 
interspectrale de signaux EEG recueillis, non seuiement ridentlfication de sources d'activites cerebrales, mais 
aussi ia localisation spatiale precise de ces sources au sein du volume cerebral. 
25 Ainsi, rinvention vise afournir un dispositif grace auquel, en particulier, la localisation de foyers ipileptlques 
peut etre realisee de lapon non traumatisante, ce dispositif n*6tant pas limite k cette application mais devant 
permettre, plus generalement, la localisation de toute source d*activite c^r^brale pathologique ou non. 

Conform§ment k I'lnventlon, le dispositif de st^rtotocallsation comports : 

(a) un ensemble d*6fectrodes destinies k 6tre dispos6es en des emplacements prSdStermines du 
S£7 scalp pour recueillir des signaux EEG, 

(b) des moyens convertisseurs pour convertir les signaux EEG recueillis sous forme numerique, et 

(c) une unite de traitement comprenant : 

- des moyens d'estimation d'une matrice interspectrale des signaux EEG recueillis, 

- des moyens de memorisation dans lesquels sont enreglstrees des informations relatives a un modete de 
35 propagation et permettant de caracteriser au molns un vecteur source dSflnissant dans le domaine 

frequentiel des relations, fonctions de parametres de position de la ou chaque source, entre un signal 
emis par une source dipolaire ou chaque source et les signaux disponibles sur le scalp, le module de 
propagation Stant un modele a trois discontinuit^s comprenant le cerveau, le crane et le scalp, et 

- des moyens de traitement de la matrice interspectrale estimee pour, d'une part, identifier les sources 
40 d'activit§s cerebrales contribuant aux signaux EEG recueillis et determiner leur nombre et, d'autre part. 

estimer les coordonnees de chaque source identifiee. 
De preference, les moyens de traitement de la matrice interspectrale comprennent des moyens 
d'estimation des vecteurs propres correspondant aux vaieurs propres de la matrice non llees aux sources 
identifiees, et des moyens de determination des coordonnees de position de sources pour lesquelles le 
45 vecteur source correspondant est orthogonal auxdits vecteurs propres. 

Le traitement de la matrice interspectrale pour ridentlfication et la localisation des sources est effeotu6 sur 
une matrice debarrassee du bruit. 

Une matrice interspectrale estimee non bruttSe peut §tre obtenu par la technique du moyennage, dans ies 
cas oCi I'activite cerebrate est provoqu6e par des stimulations, celles-ci etant repetees pour faire emerger le 
50 signal par rapport au bruit. 

Dans d'autres cas (notamment pour des sources d'activites cerebrales pathologlques), la technique du 
moyennage n'est pas toujours utilisabte. Comme indique plus loin de fapon detalllee, on peut parvenir a 
estimer une matrice interspectrale d6barrass§e de son bruit a partir de la matrice de coherence assocl6e 4 la 
matrice interspectrale estimee. 
55 L'invention sera mieux comprise a la lecture de la description falte ci-apr^s, a titre indlcatif mais non limttatif , 
en reference aux dessins annexes sur lesquels : 

- la figure 1 illustre trds sch6matiquement un mode de realisation d'un dispositif de stereolocalisation 
conforme a rinvention, 

- la figure 2 est un organigramme montrant de tapon generate les phases successives mises en oeuvre 
60 par le dispositif conforme k rinvention, 

- la figure 3 illustre un module de propagation d'ondes dans un milieu composite avec trois 
discontinuites de conductivite, et 

- ies figures 4A a 4E illustrent la localisation d'une source par mise en oeuvre de la methode du 
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projecteur orthogonal sur les donnees foumies par ia matrfce Interspectrate dSbarrassee du bruit, au 
moyen du dispositrf conforme a rinventlon. 
Dans \e dlspositif illustre par la figure 1, les signaux EEG sont recuelllis au moyen d'electrodes de scalp 10. 
puis amplifies au moyen d'amplificateurs 12, multiplexes par un multlplexeur 13 et convertis sous forme 
numerique au moyen d*un convertlsseur anaiogique-numerique 14. 5 

Avantageusement, les ampiifioateurs 12 sont du type de celui decrit dans la demande de brevet francais No 
2 564 657. 

Les signaux numerlses sont traites au moyen d'une unite de traltement 20 qui, compte-tenu du volume de 
traitements a effectuer, est de preference du type a architecture multiprocesseur organisee autour d*un bus 
commun 30. L'unlte de traltement 20 comprend : .10 

- un processGur maitre 21 assurant le sequencement des operations et le calcul final de la matrice 
interspectrale des signaux EEG recueillls, 

- un transfomriateur de Fourier rapide 22, 

- un processeur 23 effectuant des calculs de periodogrammes et interperlodogrammes, des pr^moyennages 

et des transformations logarithmiques directes et inverses, intervenant dans le calcul de la mectrice 15 
interspectrate. 

- un processeur 24 effectuant des operations de lissage optimal, 

- un processeur 25 de calcul de la matrice interspectrale debarrassee du bruit, 

- un processeur 26 de calcul de projecteur orthogonal associe a des vecteurs sources, 

- un processeur de localisation 27 effectuant le calcul de coordonnees de sources d'aotivttes cerebrales 20 
identifiees, 

- des memolres mortes (ROM) 28a destinees a ia memorisation des programmes et d'informations 
representatives de vecteurs sources pour un meddle de propagation pr6detennine entre source et 
electrodes, 

- des memolres vives (RAM) 28b pour ia memorisation temporaire de donnees, et 25 

- un circuit d 'entrees/sorties 29. 

Differents periplieriques sont connect^s a {'unite de traltement 20 dont, notamment, une imprimante 16 et 
une console de visualisation 18. 

Les phases prlncipales du programme de localisation de sources d'activites cerebrales mis en oeuvre par ie 
dispositif d6crit ci-avant sont indiqu6es sur la figure 2. qq 

Estimation de la matrice interspectrale 

La phase 101 d'estlmation de la matrice interspectrale des signaux EEG recuelllis est mise en oeuvre de la 
fapon decrlte dans la demande de brevet franpais No 87 07 600 deja crtee. 

Essentleliement, Testlmation de la matrice interspectrale est realisee par mise en oeuvre d'uh algorithme 35 
It6ratif de blanchlment et d6corr6latlon sur les transformSes de Fourier des signaux EEG recuelllis. Le 
blancliiment d'un signal est effectu^ apr^ estimation de son autospectre par la methode de lissage du 
Log-perlodogramme apr^s determination d'une fen§tre de ponderation optimale. Le passage au Log-periodo- 
gramme est de preference precede par un premoyennage du periodogramme. La decorrelatlon de deux 
signaux est effectuee apres estimation de leur coherence par la methode de lissage de rinter-periodogramme 40 
des deux signaux blanchis apres determination d'une fenetre de ponderation optimale. 

De par ses proprietes, et pour un modele de propagation donne entre sources d'aotivites cerebrales et 
electrodes disposees sur le scalp, la matrice interspectrale estimee et debarrassee de son bruit permet de 
determiner le nombre de sources decorrelees et leuf loceUisation dans le volume cerebral par rapport a un 
systeme de coordonnees. 45 

Les signaux EEG recuelllis ont un tres falbie rapport signal/bruit. Une utilisation de la matrice Interspectrale 
estimee pour identifier et iocallser des sources de fapon sQre n'est possible que dans ia mesure oO la matrice 
estimee est non bruitee. 

L'elimination du bruit peut r§sulter de la mise en oeuvre de fa technique du moyennage, c'est-a-dire en 
effectuant la moyenne de resultats obtenus par suite de stimulations repetees des sources. 50 

Lorsque la technique du moyennage n'est pas applicable, par exemple lorsque ractivitite des sources n'est 
pas liee a des stimulations exterieures controlees et n'a pas le caractere rep^titlf necessaire pour fairs 
emerger ie signal par rapport au bruit par moyennage, II est necessaire de debarrasser la matrice 
interspectrale estimee du bruit (phase 102). 

55 

Estimation de la matrice interspectrale debarrassee du bruit 

Dans rhypothese de bruits decorreles de puissance quelconque sur les electrodes, hypothese verifiee dans 
la majority des applications neurophysioioglques, la matrice Interspectrale (f; peut se mettre sous ia forme de 
la somme d'une matrice signal de rang r© inferieur a K et d'une matrice'^de bruit £ diagonaie mals non 
scalaire : = ir, * * so 

To etant la dimension du sous-espace sources (nombres de sources decorrelees) et K le nombre d'electrodes, 
tandis que 

- la notation ° signifie que Ton se trouve dans ie domaine frequentiel, 

- la notation x se rapporte aux signaux EEG recueillis, 65 
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- ta notation y se rapporte a la composante des signaux EEG provenant de la ou des sources a identifier et 
iocaliser, et 

- la notation e se rapporte k la composante de bruit des signaux EEG. ^ 
L'estimation de la matrice interspectrale dSbarassee du bruit consiste a reciiercher une estim§e r de 

5 a partir de I'estim^e ^ obtenue precedemment dans la phase 101 (la notation ^ signifie qu'll s'agit d*une 
estimSe). 

On peut montrer qu'ii suffit. pour determiner une matrice d'ordre K de fapon unique It partir de ses elements 
non diagonaux, que le rang Tq de cette matrice soit inferieur ou egal a {K/2 - 1 ) et que tout mineur (r©, ro) solt de 
rang r©. 

10 Mais, en pratique, en raison du bruit et du biais introduit par I'estimation de la matrice interspectrale, on ne 
peut affirmer quMl existe une matrice de rang r© dont ies elements non diagonaux sont identiques a ceux de 
^. De la sorte, et bien qu'en principe seuls les elements diagonaux de ^ sont affectes par le bruit, on ne 
peut recherciier qu'une estimee de a partir des elements non diagonaux de . 

Une approche par minimisation d*un critere quadratique portant sur la difference entre les elements nop 
15 diagonaux de ^ et ^ montre que I'on peut se ramener k un traitement de la matrice de coherence 
associ^e a c'est-^-<lire k la matrice dont les Sl^ments 




sont donnas par : 



25 




30 



les 



35 




40 



§tant les elements de s^, i etant rindlce de ligne et j Tindice de coionne. 
En effet, un critere quadratique de la forme 



45 




dans lequei w(i, j) est la ponderation, c'est-a-dire : 



50 




55 



peut se ramener. en remplagam -jonderation par la valeur approch^e 
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Q°P= Q /X 




Qop representant la grandeur Q normalisee de mani^re afaire apparaitre une erreur relative, d'oCi introduction 
de notion de coherence. ^ 

Le traitement mis en oeuvre sur la matrice de coherence estlmee revient done k rechercher^ dans un 
ensemble {c (r)} de matrices definies non negatives de rang r Inferieur ou egal a ro, la matrice Cm (r) qui 
minimise le critere 



Cette recherclie est effectuee a I'aide d'un algoritlime de type gradient matriciel a pas de convergence 
optimal. 

D'une fapon generate, le processus consiste a rechercher, par valeur croissante de r, dans la suite des 
{Qm (r)}, c'est- a-dire les minima de (r) , le rang r© correspondant a la plus falble valeur de r pour laquelle Q°p 
(r) est voisin d'une borne QJf (r) deflnie par : 



Le etant la largeur equivaiente de la fenetre de lissage fournie par la methode optimale d*estlmation spectrale 
des elements 



de la matrice Cx. comme decrit dans la demande de brevet franpais No 87 07 600 deja citee. Pour un falble 
rapport signal sur bruit, la borne (r) tend vers une valeur Qo^ (r) deflnie par : 
(r) = [(k-1-r)/(k-1)] QSf (r). 
SI C (ro) est ia matrice qui correspond au minimum QS" (ro) , la matrice <r debarrassee du bruit sera estlmee 
par la matrice ^ d'§l§ments : 





o 





L'algorlthme de type gradient matriciel pour la recherche du minimum sera maintenant decrit. 
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La cTifffirentielle dQ»P (r) du oritere Q"p (r), qui est nulle au minimum de ce critere. s'exprime par 



(r) 



A ^ 



A 

Ic 



10 



ce qui, compte-tenu du caractere hermitien des matrices C (r) et i, donne 



15 



. o 



20 



Si I'on designe par H (r) la matrice definle par les elements 

25 

30 

pour 1 different de j, et en raisonnant sur le numerateur Q (r) de Q^p (r) pour alleger les notations, on obtlent : 
dQ(r) = 2 tr { H(r)^dC(r)}. 

ce qui, la matrice c (r) etant factorisee sous la forme c{r) =M(r) M+(r) , conduit k la relation : 
35 dQ(r) = 2tr { H(r) M(r) dM(r)+} + 2tr { H(r) dM(r) M(rn 

la notation "tr" signifiant "trace de". 

La differentlelle dQ(r), et, par suite la differentlelle dQ°P (r) sont nulles si H(r)M(r) « 0. Pour trouver le 

mjnimum ainsi de|ini, on forme un algorithme de type gradient tel que 

5M(r) = - II H(r) M(r) avec 0 < p, 1. 
40 Pour |i suffisamment petit, raccrolssement 8Q(r) de Qfr) peut etre asslmile a sa differentielle, ce qui donne : 

5Q(r) = - 4 jitr J H fJi{r) IW(rr } . 

La trace de {HM(r) M (r)*H*} repn§sente la somme des carr6s des modules des 6l6ments de la matrice hm (r) . 

Elle est done positive ou nulle et ne peut etre nulle que si tous les elements de la matrice HM(r) sont nuls, 

c'est-a-dire pour un extrSmum de Q(r) ou, de faqon equivalente, de Q°P(r). 
45 ^ Pour une suite de^ matrices {MR(r)}, I'algorithme peut done ^tre deflnl par : 

Mh:+i{r) = Mk(r) -^k Hk Mk(r) . 
Si Ton considere la suite correspondante des valeurs Qk(r) du crlt^re, pour \im suffisamment petit, Ton 

obtient : 

Qk+i(r) - Qk(r) = - 4 |Lik tr{ Hk(r)Mk(r) Mk(rr Hk(r)*}. 
50 Cette variation est negative tant que la matrice Hk(r)Mk(r) ne se redult pas a la matrice nulle. Les valeurs {Qk 

(r)} du critere ferment done und suite decroissante qui converge vers un minimum local de Q(r)ou Qop(r). 
La Vitesse de convergence vers le minimum peut etre acceleree en recherchant le pas de convergence 

optimal, c*est-a-dire, dans la direction du gradient, le {Hk qui donne le minimum du critdre. Le p.k est determine 

de fa<?on approciiee en faisant Thypothese d'une variation localement quadratfque. A partir de ia relation 
55 ci-dessus donnant Mk+i(r), on peut ecrire : 



Mfc^l ^ktl = ^^k - l^k^ + H fc) + H ^ H+ 



En n§gligeant le terme en \Lk relatif au carre de I'erreur, on a : 

65 
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avec 



C = 



10 



< 



^CCHkii(ji^)r)=2tr{[Hk 



pour i j . 

k- l^k ^"k D^, 



O 

D 



15 



20 



De la derniire relation, en annulant la derives, II vient 



25 



ce qui permet de calculer le ^ik optima! a chaque iteration de k. 

On decrira maintenant un processus permettant dinitialiser Talgorithme afin d'eviter une eventuelle 
convergence vers une valeur erronnee par suite d'un mauvais choix de conditions initiales (c'est-a-dire ies 
ctj(r), Ce processus consiste en une methods simplifiee sejvant a estimer, par valeur inferieure, I'ordre de 
grandeur du rang recherche et a fournir une estlmee de M(r) pour initialiser ralgorithme precedent. 

La matrice c etant definie non negative peut se decomposer de fagon unique sous ie forme d'un prodult de 
matrices triangulaires avec elements positifs sur ia diagonale. Dans le cas d'une matrice de rang r inf^neur a K, 
cette decomposition prend la forme suivante : 







o o o 






o o o 

P3= BjB* 




o o o 
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Les elements non dtagonaux de c 6tant connus (identlque a ceux de c,). il suffit de determiner ies r 65 
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premiers elements diagonaux pour que la matrice soit complitement deflnie. 

En effet, le bloc diagonal Pi etant^alprs suppose connu, on peut, par un algorlthrne classigue, le 
decomposer en un produit tel que Pi = BiB"*"!. Ainsi, le blanc Bz est determine par la relation B2 P3[B i l"*" , 
ce qui permet de caiculer les elements diagonaux manquants de £ a partlr du produit matriciei : 
5 P4 = B2I . 

La determination du bloc diagonal et du rang de £ debarrassee du bruit, done la determination des r 
premiers elements diagonaux est realisee comme suit. 

On designe par { t (j)} les vecteurs colonnes de et par { les vecteurs formes par les (K-r - 1) 
dernieres composantes des (j)} , et par I la matrice formee par r vecteurs : 
10 {iHh t(2). (r)} 

En general, le vecteur ^ (r + 1) ne s'exprimera pas comme combinaison lineaire de ces r vecteurs du fait 
que la matrice \ est biaisee et brultde. On d6tinit alors 



15 ^ ^ ^ 



20 oCi les % (J) sent determines en minlmisant le crit6re 



^ o 2 

sp . . 1 w(r-fl)-> w(r-H)l 

1 w(r+l)l 



30, 

Si la matrice B est^de rang r, le vecteur param§tre correspondent au minimum est donne par 
c (r) = [ B-^ B]-^ B-^ i (rH-1j . 

En consid6rant la matrice B de rang r, les paramitres { s (j)} sont les estim6es des coefficients du (r + 1 )^^^ 
vecteur colonne de la matrice de coherence "debarrassee" du bruit enj'onction des r premiers. Pour les 
3S indices de llgne {h} oCi h est Inferieur ou egal a r, les 6I6ments diagonaux 'c hh sont done estlmes a partlr de la 
relation : 



oCi Ton supposera que % (h) est different de z^ro. 
45 Pour une premiere estimation de I'ordre, on peut ralsonner sur la valeur du crltere Q^(r). Pour une matrice 
non perturbee Q^p(r) s'annulerait pour r sup6rieur ou egal a ro. En presence d'erreur statlstique, et comme, 
plus liaut, pour le critere Q°P(r), peut determiner, en I'absence de biais, un seuil superieur Offtr) pour 
Fesp^rance du majorant du minimum du crltere Q^f'(r). Ce seuil superieur est donne par : 



55 I'Sfr+DI 



A faibie rapport signal sur bruit, cette borne tend vers 
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Dans la mesure ou Q?*(r) esttoujours un majorant du seuil d'arret, le rang estime a partir de ce critSre (qui 
ne fait intervenir qu'une sous-matrice de <f ^) correspond bien a une estimation par defaut. 

Le bloc Pi de la matrice ^tant determine comme indique au sous-paragraphe precedent et les matrices 10 
Bj et B2j,dont Tensemble definit la matrice M(r), etant calculees, on remarque que les matrices et C(r) 
(C(r) = M(r)MH-(r)) ne different que par le bloc P4 (P4 = B2 B J ). 

La matrice c(r) 6tant d^terminee entl^rement A partir des (r -h 1) premieres colonnes de c. , les lignes et 
les colonnes sont ordonnees de fapon a grouper dans ces colonnes les elements pour lesquels le coefficient 
de coherence est le plus grand. 15 

Pour ce classement 11 suffit. par example, de s'interesser a la sous-matrice Inferieure et de : 
- rechercher dans cette sousHDatrice le coefficient de coherence maximum 



tr{CJ 



1=1 



20 



et, par une permutation de Tordre des indices de la matrice, de le positionner sur la derniere ligne de la 25 
premiere colonne, 

- combiner la premiere colonne aux suivantes de fagon a faire apparaTtre des zeros sur la derniere ligne (a 
I'exception de Telement de la premiere colonne), 

- rechercher dans la matrice ainsi obtenu, privee de sa derniere ligne et de sa premiere colonne, I'element 
maximum et, par une nouveile permutation de Tordre des indices de la matrice, positionner cet element sur 30 
Tavant derniere ligne de la deuxieme colonne, et aInsi de suite. 

Ce classement conduit a une matrice triangulaire. Les Elements diagonaux de cette matrice sont ranges, de 
la gauche vers la droite, par vateur decrolssapte. Le passage par la recherche de combinajson entre les 
elements de la matrice, dans ce classement, se justifie par le fait qu'll est n6cessaire de s'assurer que la 
grandeur des coefficients de coherence ne resutte pas de la proportionnalite de deux colonnes (ou de deux 36 
lignes) de la matrice c^. 

Le processus d'estimation de la matrice interspectraie debarrassee du bruit decrit ci-avant suppose 
seulement que les bruits sont decorreles. Si Ton se place de plus dans [*hypothese bruit blanc, I'estimation du 
sous-espace source peut encore etre amelioree. 

Dans le cas ou tous les elements de la matrice sont utilises, le critere a minimiser s'ecrit alors : 40 

K 



ou C(r) est une matrice de coherence definie non negative de rang inferieur ou ega! a r et Cx la matrice de 
coherence estimee et bruitee deduite de la matrice de covariance estimee et bruitee par la relation : 



A 



y /a a 



Sous I'hypothese bruit blanc, la matrice C(r) peut s'e 
C(r) = M(r) M+(r) + oHr) D i 



so 
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ou M(r) M+(r) represents la matrice de coherence debarrass^e du bruit au rang r, a^(r) la puissance du bruit 
blanc au rang r et Dx la matrice diagonale definle par : 



10 25^ = 



15 



0 



20 



25 



Dans ces conditions, il est possible de determiner le aHv) qui minimise Q'«'P(r), 
tr tC^Cr) -M(r)M"*"(r)} D"^ 



cr(r) = 



A. o 

tr {Df } 



Ainsi, si I'on pose 
30 c , = Cx-oHr) Dx 

le rang ro optimal est d§fini en minimlsant le critere : 
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au iieu du critere Q°P(r). Le reste du processus decrit preti-^demment reste Int^hange et s'applique 
directement. ^ 

Le calcul iteratif sur r et sur a h permet didentifler sous-esr. source en determinant To . d'estimer ia 
puissance du bruit blanc et, de plus, fournit une estimation de la ^Ttatdce de covarlance debarrassee du bruit. 

L'estimation de la matrice Interspectrale d§barrassee du bruit permet d*evaluer le nombre To de sources 
decorr^l^es (cette approctie pouvant s'6tendre au cas de sources correl§es) et de fournir, par le calcul des 
valeurs propres associees, une information relative k la puissance d'emission de ces sources. Mais 
Testimatlon de la matrice interspectrale debarrassee du bruit permet aussi la localisation des sources, dans la 
mesure ou un models de propagation entre sources et electrodes a ete definl. 

Definition d'un modeie de propagation 

Le modeie utilise ici est un modeie a trois discontinuites de conductivlte electrique schematise par trois 



spheres concentriques de rayons respectlfs ai, 
cephalo-rachldlen, le crSne et le scalp (figure 3). 

On se place dans I'hypothese d'une charge ponctt>^ 
spheres. Les dlscontinuit6s de conductivity sont tre ' 
100 entre la conductivlte du cr&ne et les conduce; 
Tune de Tautre, et puisque la conductivity 04 ^ Textene 



a2. as representant le cerveau avec le liquide 

:>lacee en un point Mo a distance b du centre O des 
. ortante, puisque Ton peut atteindre un rapport de 
r rsA et as du cerveau et du scalp qui sont voisines 
, :-u scalp est trds differente de celle du scalp.(a4 « 0) . 
Du fait de ces importante discontinuites de conductivite, la contribution en sortie des ondes transmises 
intermediaires peut etre negligee. On ne considere done, en sortie, que la contribution des ondes Dnm ayant 
subi n reflexions entre les deux couches 31, 32 du cote interne et m reflexions entre les deux couches 32, 33 
du cote externe. 

La position Mnm sur le segment OMo de Timage (source virtuelle) dont est Issue Tonde Dnm est donnee par 
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Quant II i'amplKude C^nm de I'onde Dnm, elle est donnee par 



Le potentiel resultant k la surface de ia sphere exterieure s*exprime done par : 

P etant Tangle forme par le rayon passant par le point consldere k la surface de la sphere par rapport au rayon 
OMo passant par la source, et^o (P, OMnm) etant le potentiel en surface d'une sphere du a une charge unite 
placee en Mnm a la distance OMnm du centre de la sphere pour un milieu homogene. 
(P, OMnm) s'exprime par : 



10 



26 



30 



35 



40 



45 



as = rayon de la sphere exterieure. 

La double somme infinie donnant le potentiel Yr (P. b) peut §tre slnipliflee par approximation par une 50 
integrate continue de la fornrie 



^ o 
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LMntegrale ci-dessus peut etre calculee par la methode des trapezes en discretlsant la variable t, ce qui 
revient a placer des charges regulierement espacees entre le centre O de la sphere et la position Mp de la 
charge reelle avec des poids donnes par la fonctlon de ponderation : 



Le module decrlt ci-avant pour une charge ponctuelle est facilement etendu a une source dipolaire qui offre 
une meilleure representation d'une source d'activites cerebrales. Le passage a la source dipolaire est obtenu 
nunneriquement en sommant les potentiels crees par deux charges--^ et 4- (ou 1 caracterise le moment 
dipolatre et e leur distance) placees a distance (b - e) et b du centre des spheres avec e tres petit devant b, et 
en remarquant que, dans le cas d'une sphere conductrice le potentiel a la surface de celle-ci peut s'interpreter 
comme la somme du potential de deux dipdies dans I'espace libre de m§me conductivity solt, finalement, k 
partir de la charge en Mnm : 

[^'oCP,!!) , -H'o(P,ii-e)] et [V2^o(P, -PA ), -V2To(P, )]. 

A titre d'exemple, en caiculant par la methode des trapses I'lnt^grale definle plus haut donnant une 
approximation du potentiel resultant (P, b), en utilisant cinq traptees (charges virtuelles plac6es k 3/16, 
1/8. 1/4, 1/2 et 3/4 de la distance OMo k partir de O, on obtient pour le potentiel cree par un dipdie reel : 



Le modele de propagation qui vient d'§tre d^crit permet done de d§finlr un modele de vecteur source S (cp) 
dent les composantes repr^sentent les valeurs des potentiels recuelllis (a une constante pi^s) sur les K 
electrodes en fonction des coordonnees de la source, coordonnees exprimees par un vectt>(!r de position q) 
relatif a la source. En d'autres termes, le vecteur source S{(p) caracterise la fonction de trrtnsfert dans le 
domaine frequentiel entre la posrtlon de la source et les signaux EEG recuelllis sur les electrodes. 

Le modele de vecteur source S((p) correspondant au modele de propagation selectlonne est enregistre en 
memoire ROM du dispositif de traitement pour §tre utilise auxfins de localisation des sources Identifiees dans 
la matrlce Interspectrale dSbarrassee du bruit estimee Ty. 

Localisation des sources 

De preference, la localisation des sources est effectuee par mise en oeuvre de la methode du projecteur 
orthogonal, ou encore la methode du gohlom^tre g^neraiisee (phase 103). 

Ayant constats que pour le modele dipolaire, le vecteur source^est independant de la frequence, on 
localisera la source, non pour chaque plan frequentiel de la matrice , mais pour une matrlce ry"* obtenue 
par moyennage sur les differents plans frequentlels. 

La matrice possede (K - ro) valeurs propres qui sont voisines de zero (si Testlmation etait parfalte, ces 
(K - ro) valeurs propres seralent nulles). A ces (K - ro) valeurs propres correspondent des vecteurs propres 
{ S]} qui sont orthogonaux aux vecteurs sources { Ii} qui ferment le sous-espace source de dimension r©. 

On designe par p((p) la fonction obtenue en projetant le meddle du vecteur source S(<p)sur les valeurs 
propres { -Sj}, avec normalisation par S(<p): 
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f^.f etant ie conJugu6 transpose de ^j. 

A chaque fois que Ie vecteur de position <p sera egal au vectegr de position (p z d'une source, ia fonction io 
p(cp) dolt etre nulle a cause des fluctuations statistlques. C'est done la recherclie des minima de cette fonction 
(ou les maxima de son inverse) qui donnera les vecteurs de position reciierches (pi des sources. On notera 
que, dans I'expression ci-dessus de p((p), ia normalisation par |s ((p)| est introduite pourtenir compte du fart 
que ia norme des vecteurs souces n*est pas obligatolrement constante. 

La recherche de minima de p((f>) est effectuee en balayant {'ensemble des valeurs du vecteur de position . is 
En variante, ta localisation des sources pourra etre effectuee en utilisant Ie sous-espace source. Une 
mithode peut consister alors k minimlser ia distance entre 



et etant les valeurs propres et vecteurs propres du sous-espace source, et 



25 



30 



35 



est la densite spectraie de puissance de la source i. 

Toutefois, par rapport a la methode du projecteur orthogonal, cette methode presente Ie desavantage de 
fatre porter la recherche des parametres de position sur tous les parametres a la fois. 

On notera que ia methode de localisation par projecteur orthogonal decrlte ci-avant en liaison avec le 40 
modele de propagation envisage impose une disposition des electrodes dans un meme plan et suppose une 
orientation radiale du dipole representant la source. Le choix du modele dipolaire radial est ici justifie en raison 
de rorientation radiale des cellules pyramidaies au niveau du cortex. 

La disposition des Electrodes dans un plan permet une localisation uniquement dans ce plan. 

Une localisation en trois dimensions est possible en effectuant des localisations dans plusleurs plans 45 
paralleles en liaison avec une cartographie donnant dans chaque plan une id6e sur les niveaux de puissance 
d'emission. 

La localisation en trols dimensions peut etre directement obtenue en eiaborant un modele de propagation 
plus complexe faisant intervenir un parametre de position supplementaire et une orientation quelconque du 
(ou des) dipdies. so 

Exempie 

Les figures 4A a 4E se rapportent a un exempie de mise en oeuvre du dispositif seion IMnventlon pour ia 
localisation d'un foyer epileptique. 

La figure 4A montre les emplacements des 16 electrodes utillsees. S5 

La figure 4B montre les signaux anaiogiques EEG recueillis sur les 16 Electrodes. 

Une matrice interspectraie n; a ete estimEe a partir des signaux EEG recueillis par les Electrodes 
coplanalres No 4, 10, 9, 13, 14 et 8. ^ 

La diagonaiisatioh de ta matrice de coherence associee a la matrice Interspectraie estlmee clans la 
bande delta (0,1 a 4 Hz) fait apparaitre six valeurs propres non nulles (autant que de capteurs) dont deux 60 
paraissent prepqnderantes. 

Une matrice f y debarrassee du bruit a ete estlmee a partir de comme decrit plus haut. De Testimatlon 
de iPy, 11 ressort que la matrice a un rang egal a deux. La diagonal Isation de la matrice de coherence cy 
associee a tFy fait apparaitre deux valeurs propres distinctes auxquelles peuvent etre associees deux sources 
decorrEiees. 65 
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La localisation des deux sources a ete effectuee par recherche des maxima de la fonction 1/p{q)), ce qui est 
equivalent a la recherche des minima de p(<p), en utilisant la mSthode du projecteur orthogonal avec deux 
vecteurs source elabor§s k partlr du models de propagation a trols discontinuitSs de conductivity, comme 
decrit plus haut. ^ 
5 La figure 4C montre la cartographie realises a partlr de la matrice interspectraie estimee n^, dans les 
bandes alpha lent (7,5 a 10 Hz), beta rapide (18 a 40 Hz), theta (4 a 7,5 Hz), beta lent (12,5 a 18 Hz), delta (0,1 
a 4 Hz) at alpha rapide (10 a 12,5 Hz). 

Les cartographies dans la bande theta et delta indiquent la presence d'un foyer epileptique dans fa region 
temporale drolte. 

10 La figure 4D montre la variation de la fonction 1/p(cp) dans le plan defini par les electrodes selectionnees. 
Deux maxima de la fonction 1/p(<p) apparaissent sur la figure 4D, les positions de ces maxima 6tant donnees 
sur la figure 4E par une croIx et un asterisque. Sur la figure 4E, l*axe horizontal passant par le centre O 
correspond a celui des Electrodes No 14 et 10. l-a croix correspond a la iocaiisatlon du foyer epileptique. tandis 
que rast^risque refidte vralsemblablement une "source" relative a des mouvements oculaires. 

15 



Revendications 

20 

1. Dispositif de stereolocalisation spatiale de sources d'activites cerebrates, caracterise en ce qu'il 
comporte : 

(a) un ensemble d'electrodes (10) destinees a etre disposees en des emplacements 
predetermines du scalp pour recueiilir des signaux EEG, 

25 (b) des moyens convertisseurs (14) pour convertir les signaux EEG recueltlls sous forme 

numerlque, et 

(c) une unite de traltement comprenant : 
-des moyens d'estimation optlmale d*une matrice interspectraie des signaux EEG recueillis, 
- des moyens de memorisation dans lesqueis sent enregistrees des informations relatives a un 

30 modele de propagation et permettant de caracteriser au molns un vecteur source def inissant des 

relations, fonctlons de paramdtres de position de la ou chaque source, dans le domaine frequentle!, 
entre un signal emis par une source dipolaire ou chaque source et les signaux disponlbles sur le 
scalp, le module de propagation Etant un module k trols discontinuites comprenant le cerveau, le 
crane et le scalp, et 

55 . des moyens de traltement de la matrice interspectraie estlm§e pour, d'une part, identifier les 

sources d*activites cerebrates contribuant aux signaux EEG recueillis et determiner leur nombre et, 
d*autre part, estimer les coordonnees de chaque source identifi§e. 

2. Dispositif de stereolocalisation selon la revendication 1, caracterise en ce que les moyens de 
traltement de la matrice interspectraie comprennent des moyens d'estimation, sur la matrica 

40 interspectraie debarrassee du bruit, des vecteurs propres correspondant aux valeurs propres de la 

matrice non llees aux sources identifiees, et des moyens de determination des coordonnees de position 
de source pour lesquelies le vecteur source correspondant est orthogonal auxdits vecteurs propres. 

3. Dispositif de stereolocalisation selon Tune queiconque des revendications 1 et 2, caracterise en ce 
que les moyens de d'estimation de la matrice Interspectraie dibarrassee du bruit, a partir de la matrice de 

45 coherence associee^ la matrice interspectraie. 
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